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ABSTRACT 

Context. The globular clusters of our Galaxy have been found to lie close to a plane in the logRg, logo", SB^ space (see Djorgovski 
1995), on the continuation of the Fundamental Plane that is known to characterize the global properties of early-type galaxies. Of 
course, there is no apparent reason why such physically different self-gravitating systems should follow the same scaling law. 
Aims. We reexamine the issue by focusing on a sample of 48 globular clusters selected in terms of homogeneity criteria for the 
photometric data available from the literature. 

Methods. We perform a model-independent analysis of surface brightness profiles and distance moduli, estimating error bars and 
studying selection effects with robust non-parametric statistical tests. 

Results. We determine the values of the coefficients that define the Fundamental Plane and their error bars and show that the scatter 
from the Fundamental Plane relation is likely to be intrinsic, i.e. not due to measurement errors only. Curiously, we find that in the 
standard Fundamental Plane coordinates the set of points for our sample occupies a rather slim, axisymmetric, cylindrical region of 
parameter space, which suggests that the relevant scaling relation might be around a line, rather than a plane, confirming a result noted 
earlier. This is likely to be the origin of the difficulties in the fit by a plane often mentioned in previous investigations. In addition, 
such a Fundamental Line relation would imply a pure photometric scaling law relating luminosity to the effective radius which might 
be tested on wider samples and on extra-galactic globular cluster systems. As to the residuals from the Fundamental Plane relation, 
we find a correlation of the deviations from the plane with the central slope of the surface brightness profile. No other statistically 
significant correlations are identified. Finally, given the constraint imposed by the virial theorem, we study the distribution of the 
values of the quantity Kv/(M/L) (virial coefficient divided by the relevant mass-to-light ratio); the distribution of the logarithms, 
reconstructed through kernel density estimation methods, shows evidence for bimodality, which suggests that the galactic globular 
cluster system may be composed of at least two dynamically different populations. Yet, these populations do not appear to reflect the 
standard dichotomy between disk and halo clusters. 

Key words. Galaxy: globular clusters: general - Galaxy: structure 
1 . Introduction quasars (e.g.. iHamilton et al]|2006l) . Comparisons of the FP re- 



T • „ , JT^ \ . [tt: . . „ ^ — . lation for differ ent systems are often pr esent in the literature : 

It IS wel known ( [Dressier et alJl987| ;| Diorsovski & Davis|l987|) for example, see ^urstein et alJ ([19971) and lMartini & HolJIOOl . 

^ that early-type galaxies occupy a region in the vicinity of the ^j^^ ^j^^ ^^^^^^^ 

problem of the mass gap between high- 

^ plane log/?, = aloga-„ + /3SB, +ym the three-dnnensional ^^^^ low-mass galaxies and related classification is- 

^ parameter space defined in terms of the projected effective ra- ^^^^ introduced by Bender et al. ( 1992i). 

rrt diuajKp, of the central velocity dispersion cTn, and of the mean j . , r i /-.^ . i. c j i 

, . , r,r, n 1 r- J i r,i zT-nx i -1 • In particular, for galactic GCs it has been found by 

■ " " surface brightness Such Fundamental Plane (FP) exhibits a irr: — i N inng-k >u . i-n w u- u j 

., . , , , r , m ■ , . • IDiorgovskil(ll995D that an FP relationship based on core param- 

tilt. I.e. the measured values of the coefficients that characterize ~ , ... , , , ^ u ■ \ u u j . 

, , , ^ , . . , , eters (core radii and central surface brightness) holds and that 

the plane differ sigmhcantly from the naive expectations based • • ^ . -.u i i. i i . i ■ • • i- j 

, ,• • .c 1 ••11 r„, • ^ , It IS consistent with globular cluster cores being virialized sys- 

on the application of the virial theorem. The properties of the , . i j . . ,^,t .• u^i c 

r-r^ r. , , , ■ ■ r , , • • tcms With a umvcrsal and constant M/L ratio , while for quanti- 

FP reflect partly the characteristics of the stellar populations in ,• ^ j » a ^- /n j m \ r-n j 

, I , • ■ • • .c • ties referred to the effective radius (Kf and 5B-,) the FP does not 

these systems and partly some systematic deviations from strict .1 •i.i 1. r w / 

■' , / |_ . — n i^oo^i . emerge as sharply, possibly because of error correlations (see 

structural homologv (e.g.. see Bertin et alJ2002l) . , aa-^- 11 c » mu 1 n t r» • r^r^ 

- . ^ , , , '■ 7^^ r, r^r, . also additional remarks in Sect. Blbelow). In reality, since GCs 

It is often stated that the universality of the FP extends u 1 . n . •.u u 

„ • , , , • J- , , • , ,• are non-homologous stellar systems with rather homogeneous 

well outside the domain of early-type galaxies and applies ^ „ w •^ u 1. • • •r u j . 

, , , , , , /I J- stellar populations, it would be surprising if they happened to 

to many other systems, such as globular clusters (hereafter „ .• »• f^u i-n f 1 . 1 • u- u 

rrr — — nnr^TT?! ITT T , .• ^^^^ r. , n follow a continuation of the FP of early-type galaxies, which ap- 

GCs; see PiorgovskJ 1995U McLaughlin 2000; BarmbyetalJ ^ , ■ , a -.u A-tr , 1 ■ 

u„^Ji , — — , — . , „ f„ . , , ^„„,, l pear to be associated with very different structural properties, in 

20071) . cluster spheroids (Zaritskv et al. 2006), open clusters f f, , ^ , ^. 1 „ » » 

'i _ ' — - _. — ii i: I ,,. y , , . terms of homology, stellar populations, and dark matter content. 

(^Bonatto & Bica 2005), X-ray emitting gas in elliptical galaxies t.^ . •j. 11 r , 

,T^^ 1 1 D n 1 1 r , ■ / R ■ — TTI In this papcr we thus reconsider the problems of the exis- 

(Diehl & Statler 2005), clusters of galaxies (see ILanzoni et alj , , ^ r r-r. r /-^^ r l,^t■^^ 

r,„„ . A £ li. ■ \ A — 1^ T- — 7^ tence and of the nature of the FP for the GCs of the Milky 

12004 , and references therein), and even more exotic obiects, ^ , ... , , • , ,, ,• , 

r • 111.1 / m — T- :~i ^rvrvrl A Way Galaxy. Although galactic GCs are extremely well studied 

such as supermassive black holes (e.g.. lHopkins et al. 2007 ) and , .■' „ . . j- 1 • 1 
° — ^ objects, for a proper statistical investigation of their structural 

' Some authors employ r„ R,„ or n, to refer to the projected efi^ective properties the use of the entire set of data available from the 

(half-light) radius. In our notation r, refers to the projected efi'ective Uterature suffers from several limitations and is therefore prob- 

angular radius, while to the projected effective linear radius. lematic. In fact, since the GCs are basically distributed over the 



2 



Pasquato & Berlin: Fundamental Plane of Globular Clusters 



entire sky, the relevant photometric and spectroscopic data that 
have been collected come from different telescopes and different 
instruments. In addition, GCs are located at different distances 
from the Sun and, in the sky, at very different ranges of galac- 
tic longitude and latitude, thus suffering from different amounts 
of reddening. Thus different groups of GCs suffer from different 
uncertainties in distance determination, which has a major im- 
pact on the determination of the intrinsic properties that enter the 
FP. All this makes the entire set of data highly inhomogeneous 
and thus not well suited for a proper statistical investigation. 

As a preliminary step for a satisfactory statistical investiga- 
tion, we have thus decided to identify a sample of galactic GCs 
that, in our view, has optimal characteristics, being sufficiently 
large while being associated with sufficiently homogeneous data 
in relation to their photometric properties. Our study of the FP 
has thus been restricted to such selected sample. Sample selec- 
tion is a delicate issue, so a thorough discussion of this problem 
is carried out in the following. 

The paper is organized as follows. In Sect. 2 we describe 
the criteria that have led us to identify our optimal sample of 
GCs and then summarize its properties. In Sect. 3 we proceed 
to study, for such sample, the possible existence of the FP in the 
relevant parameter space and show that, in the standard FP co- 
ordinates, our sample occupies a rather slim cylindrical region, 
which suggests that the scaling relation might be around a line, 
rather than a plane, a result alre ady noted in previou s investiga- 
tions (Bellazzini 1998, see also iBurstein et alJ[T997h . In Sect. 4 
we discuss two unexpected features, i.e. a residual correlation 
with a parameter characterizing the central slope of the bright- 
ness profile and a bimodality in the distribution of the effective 
virial coefficients Kyi (MIL), reconstructed by means of kernel 
density estimation methods. In Sect. 5 we discuss the selection 
effects on our sample in comparison with those present in pre- 
vious investigations (in particular, in Diorgovski 1995). Finally, 
in Sect. 6 we summarize the results obtained in the paper and 
present discussion and conclusions. 



2. The adopted sample 

Since our primary objective is to test whether indeed the galactic 
GCs follow a continuation of the FP of early-type galaxies, we 
have decided to refer to the standard global photometric param- 
eters, i.e. to the effective radius and to the related mean surface 
brightness (Rg and SB^), which are commonly used in the study 
of the FP of early-type galaxies (for a complete definition, see 
Bertin et al. 2002) and to ignore other options (such as the use of 
core parameters, /?() and SB()\ see Djorgovski 1995). The criteria 
for the identification of our sample, extracted from the data set 
available from the literature, are largely dictated by the require- 
ment of accurate distance determination. 

In the following four subsections, we address separately the 
main ingredients required for an accurate determination of the 
parameters entering in the FP relation, i.e. distance moduli, red- 
dening, surface brightness profiles, and line-of-sight velocity 
dispersions. We will then complete the Section by identifying 
our sample. 

2.1. Distance moduli 

The two most recent papers that contain an extensive set 
of galactic globular clusters for which distance moduli have 
been determined with a uniform method (using the Zero Age 
Horizontal Branch (ZAHB) as a standard candle) and have been 



assigned individual error bars are the article by iFerraro et al.l 
(1999"), who studied and reduced uniformly colour-magnitude 
diagrams com ing from several sources , for a set of 61 GCs 
and that by iRecio-Blanco et akl (l2005h . based on an exten- 
sive photometric data set from HST for 72 GCs (iPiotto et al.l 
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120021) ■ The samples overlap partially. Different procedures are 
used in the two papers to deal with data reduction, partic- 
ularly with respect to the determination of the ZAHB level. 
Recio-Blanco et al. (2005) compare their distance moduli whith 
those of iFerraro et al.. (.1999) on the intersection of the two sam- 
ples, a nd find them on average larger than those of lFerraro et al.l 
?1999) with a mean offset of 0.09 mag. Moreover, an rms scatter 
of 0.17 mag is found for the differences, with some values being 
as discrepant as 0.4 mag. Therefore, in this paper, we decided 
to consider the properties of a combined sample, which contains 
48 GCs (see Sect. 12. 5t , but also to check the behaviour of two 
sub-samples studied separately (se e Appendix). 

In their Table 2, IFerraro et af] ([1999,) fist ZAHB magnitude 
levels, with the related photometric error bars, and two sets of 
true distance moduli (\m - MIo), based on two different assump- 
tions on the enhanced abundance of a-elements (the so-called 
[a/Fe] enhancement) of GC stars with respect to solar abun- 
dances. Column 7 of that Table lists the distance moduli ob- 
tained by assuming solar abu ndances (i.e., with no [a/F e] en- 
hancement) and [Fe/H] from Carret ta & GrattonI (Il997h . while 
column 8 lists the distanc e moduli obtained by adopting \a/Fe] 
values from the literature dSalaris & Cassisi ' 1996; C arnevll 19961) 
and by mimicking ff-enhanced isochrones by stand ard scaled so- 
lar models of suitable metallicity, as proposed by ISalaris et al.l 
(il993ft . The authors state that distance moduli "are affected by 
many uncertainties (namely, the evaluation of the ZAHB level, 
the zero point and dependence on metallicity of the ZAHB level, 
reddening, etc.)" so that "the global uncertainty affecting the dis- 
tance moduli listed in Table 2 cannot be less than 0.2 mag". 

If we linearly interpolate the distance moduli as a function 
of [a/Fe] on the interval between and the value of [a/Fe] 
adopted for the construction of column 8, we can take into ac- 
count the uncertainty on [a/Fe] by adding to the distance mod- 
uli error the product between the relative error in [a/Fe] and 
the difference between the distance moduli listed in Table 8 and 
those in Table 7. This appears to be reasonable. In view of the 
sparseness and hete rogeneity of [a/Fe] data, as noted also by 
Ferraro e t alj (Il999l) . we estimate the relative error on this quan- 
tity to be approximately 50%. With this device, we can take 
into account this factor in the uncertainty on dista nce moduh. 
Keepin g in mind the level of 0.2 mag suggested by Ferraro et al.l 
( Il999h . we decided to proceed by linearly transforming the er- 
ror bars obtained above, so that the first quartile of the error 
bar size distribution was rescaled into 0.2 mag and the third 
quartile into 0.3 mag. This procedure is both reasonable and 
statistically robust, because it is not affected by outliers in the 
error bar size distribution and turns out to be compatible with 
the scatter and systematic differences that arise when compar- 
ing distances obtained from different standard candles, such as 
the Main Sequence or the White Dwarf Cooling Sequence (for 
a distance measuremen t of NGC 104 based on the White Dwarf 
Cooling Sequence, see IZoccali et al. 19991) , or ev en when com- 
paring the distance r nodul i of lFerraro et al.l (119991) with those of 
iRecio-Blanco et alj (l2005h . 

The distance moduli o fH ecio-Blanco et al. I (l2005h are prob- 
ably of higher quality with respect to those of 'Ferrar o et al.l 
( 1999), because they are based on uniform HST photometry, at 
least as far as random error bars are concerned. For the inter- 
section of the two samples our adopted uncertainties are there- 
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fore obtained by summing in quadrature to the error estimates by 
iReci o-Blanc o et all (|2005) the semi-difFerence between the dis- 
ta nce moduli obtained by Recio-Blanco et al. (2005) and those 
bv lFerraro et a l. ( 1999) . For the remaining clusters with distance 
moduli only in .Recio-Blanco et alj (120051) . we assume a constant 
difference equal to the mean difference quoted above, and use it 
to do an analogous quadrature sum to get the adopted error bars. 

Our conservative procedure may tend to overestimate the er- 
ror bars on distance moduli. Still we will see that the formal 
reduced for the FP fit will turn out to be larger than unity. 

2.2. Reddening 

iRecio-Blanco et all (|2005[ ) provide independent measurements 
of the reddening affecting the GCs in their sample, which we av- 
erage d with those reported in the McMaste r Catalogue (iHarris 
1 19961). For the GCs not considered by iRecio-Blanco et aL 
(2005), we just adopted the values reported in the latter 
Catalogue. On the values of reddening, we assigned error bars 
as suggested by Harris: "The typical uncertainty in the redden- 
ing for any cluster is on the order of 10 percent, i.e. 6[E(B - 
V)] = Q.IE(B - V)". Note that the appa rent distance moduli 
determined by iRecio-Blanco et al.1 (|2005) refer to the F555W 
HST band; to de-redden them, we thus applied the relevant 
reddening/extinction r elationship, obtained from Table 12 of 
iHoltzmanetai] (119951) through a linear fit of the foimAfsssw - 
kE{B — V) and differs from the usual Ay - 3.\E{B - V) essen- 
tially in the value of the proportionality coefficient k. 

2.3. Surface brightness profiles: apparent magnitudes and 
effective radii 

To calculate apparent effective radii and apparent V- 
band magnitudes (or equivalently, the values of mean V- 
band surface brightness within the effective radius, SB,,) we 
started fro m the circul a rized surface brightness profiles re- 
ported b y 'Trager et al.' ('1995V A more recent compilation 
dMcLau ghlin & van der Marel 2005) of GC surface brightness 
profiles and structural parameters is available, comprising both 
Milky Way and extragalactic GCs, but, for the galactic GCs, 
it is based on fitting various models to the same set of data 
by 'Trag er et alj (|j_995'). Our model independent approach thus 
prompted us t o pr oceed independently, basing our analysis di- 
rectly on the Trager et alj d 19951) profiles. These profiles are 
based on data from different sources and do not quote explicit er- 
ror estimates. On the other hand, the authors provide some qual- 
ity rating of the data, based on the reliability of the data-set from 
which they were extracted. These weights, w, ranging from 0.0 
(totally unreliable data-set) to 1.0 (high quality data-set) were 
admittedly "assigned by eye" (Trager, private communication). 
Relative error bars might be introduced as cr oc 1 / -\/w, but this 
would still be somewhat arbitrary. Eventually, to set the magni- 
tude of the error bar for each surface brightness data-point, we 
resorted to the (halved) residua ls to the Chebyshe v polynomial 
fits to the profiles provided by ITrager et al.l d 19951) . because we 
wished to avoid relying on a subjective assessment of data qual- 
ity. This choice has its own limitations, because such residuals 
depend on the particular form of the fitting function, are partly 
correlated, and may include as an error term some spurious de- 
viations (e.g., deviations associated with the presence of a bright 
star). 

The physical parameters that we wish to derive require inte- 
gration of the profiles on the radial coordinate. Since the points 



are not evenly spaced in radius and are rather noisy, we have 
decided to smooth and interpolate the data by cubic splines, 
prior to integration (see also Trager etal. 1995). We then pro- 
ceeded to integrate using different methods (e.g., trapezoidal and 
Simpson's rule) and checked that the results were stable. The 
integration error was obtained by comparing results derived by 
changing the number of interpolating points and controlled to 
be at least one order of magnitude smaller than other sources of 
error 

The surface brightness profiles were extrapolated beyond the 
last available data-point, which is set by the merging of the pro- 
file into the background sky brightness and so differs from one 
GC to another To extrapolate, we used a linear relation between 
surface brightness and angular radius. We tested a number of al- 
ternatives for an optimal extrapolation and noted that the result 
did not change the derived total magnitude by more than 0.05 
mag. 

The apparent effective radii were then determined in a 
straightforward manner by building a curve of growth through 
integration of the profile and finding the radius at which the inte- 
grated luminosity equals half of the total luminosity. Obviously, 
since uncertainty on the total magnitude propagates to the effec- 
tive radius, this factor has the potential of introducing an anti- 
correlation between uncertai nties on the two quantities (see also 
Sect. 3.1 o f iTreu et al.llT999h . The importance of this error term, 
which in the error budget for the effective radius competes with 
the error deriving from the uncertainty on the parameters of the 
interpolation curve, depends strongly on the steepness of the 
growth curve in the vicinity of the effective radius. 

Fig. [T] shows a comparison of the values of effective radii 
obtained in this paper with those recorded in the M cMaster 
Catalogue (Harris 1996). Following iTrager et al.l (1 19951) . we re- 
placed points characterized by excessive dis crepancies wi th val- 
ues from the literature (in our case from lHarrislll996h . This 
was performed by replacing only those data that we found to 
be outliers in the distribution of differences. For this purpose, 
we defined as outher a point which lies outside the interval 
[Qi - (63 - 20,23 + (G3 - Qi)l with Qi and Q3 being the 
first and third quartile of the data distribution; this is slightly 
more conservative with respect to the more common definition 
based on [gi - 1.5(fi3 - gi), Q3 + 1-5(23 - Qi)]- 

Our angular radii are slightly larger than those listed in the 
McMaster Catalogue, showing a median difference of 0.04 in 
the logarithm. The standard deviation of the differences is some- 
what larger, at 0.10, but still compatible with our estimated error 
bars on effect ive radii. Similarly, our magnitudes exceed those of 
lHarrisl(ll996h bv 0.07 mag in the median (i.e., we find clusters to 
be systematically dimmer), and the differences show a standard 
deviation of 0.21 mag. This scatter is slightly higher than that 
expected based on error estimates for our sample alone. 

Unfortunately, the surface brightness profiles presented by 
ITrager et al.l (1 19951) do not allow us to carry out the analysis in 
different colors. Therefore, we will be unable to address the pos- 
sible impact of mass segregation on the FP scatter. 

2.4. Velocity dispersions 

We took the line-of-sight velocity dispersions from 
iPrvor & MevlanI (Il993h . which are complete with uncertainties 
estimated by the authors. The term central velocity dispersion 
(and the corresponding subscript 0) would be misleading in this 
context; so, while we have used it when referring to the FP of 
elliptical galaxies, we will drop it in the following, because the 
data set contains some dispersions that are not strictly referred 



4 



Pasquato & Berlin: Fundamental Plane of Globular Clusters 



2.5. The sample 



9- ° 
ra o - 

Cl c\j 




50 



100 150 
re(arcsec, Harris 1 996) 



200 



250 



Fig. 1. Comparison of apparent effective radii reconstructed 
in this paper with those reported in the McMaster Catalogue 
(lHarrislll996.) . A 1 : 1 line has been drawn to guide the eye. 



to the center of the cluster We did not try to rescale these mea- 
surements to a common aperture size, because this woul d be, 
tosome extent, model de pendent (e.g., if we made use of IKind 
(119661) models; see also iMcLaughlin & van der Marel (l2005h 
and Table 13 therein for observed central velocity dispersions 
corrected by using different models). In the adopted data set 
two kinds of velocity dispersion measurements are mixed, 
from integrated light spectroscopy and from the kinematics of 
individual stars. On a total of 56 GCs, 31 have integrated light 
measurements and 35 have stellar kinematic measurements, 
while only 10 have both. 

Integrated light measurements make use of a small (several 
arcseconds) spectroscopic slit centered on the GC. In contrast, 
samples of stars used for stellar kinematics often span a much 
wider angle (up to several arcminutes). In addition, integrated 
light spectroscopy cannot discriminate between the cluster and 
background or foreground sources, nor can it exclude fast mov- 
ing binary stars which may exaggerate the observed velocity dis- 
persion. Moreover, the light collected by the spectroscopic slit 
may be dominated by a small number of bright giants, therefore 
amounting to poor statistics. On the other hand, stellar samples 
used for kinematics are usually selected carefully for member- 
ship and tested against binaries through repeated observations, 
but may be rather small and unevenly distributed in the cluster; 
they also are often biased towards highly luminous stellar types 
(mainly red giants). All these factors contribute to producing a 
systematic difference between these two kinds of measurements, 
which we have tested on the subsample of clusters having both 
kinds of measurement and found to be significant. In Sect. 14.21 
devoted to the analysis of the distribution of virial coefficients, 
we will restrict our study to the relatively homogeneous sub- 
sample of GCs for which the velocity dispersion is given by 
integrated light measurements, because the virial coefficient is 
proportional to the square of the velocity dispersion and sample 
homogeneity is crucial. For the rest of the paper, we repeated our 
calculations on the restricted subsample and checked that the re- 
sults are in general agreement with those obtained on the larger, 
unrestricted sample. 



The McMa ster Catalog ue of Milky Way Globular Cluster 
Parameters (' Harrislll996 ') lists 151 globular clusters. The sam- 
ple we selected is much smaller, comprising only 48 globular 
clusters. Moreover, as explained below, data inhomogeneity has 
forced us to consider separately also two smaller subsamples for 
which homogeneous measurements of distance moduli are avail- 
able. 

Since we are aiming at a model-independent approach, we 
have carefully excluded from our sources those studies that pro- 
vide exclusively model-dependent quantities (in particular, sur- 
face brightness profiles from King model fits). Since we wish to 
determine error bars on the FP coefficients, we have restricted 
our attention to data with objective estimates of uncertainties. 
Finally, we have given priority to homogeneous data-sets (i.e. 
data from the same authors, taken with the same instruments, 
and reduced in a similar fashion). 

In conclusion, our sample has been selected as the in- 
tersection of three data-sets: one with well-deter mined dis- 

(119991) 



tance moduli obtai ned by 
iRecio-Blanco et all (12005^ 



mergmg iFerraro et al.1 (119991) and 
one with well-determined line-of- 



sight velocity dispersions dPrvor & Meylan 1993), and one with 
surface brightness profiles (Trag er et al.lll995l) from which we 
could recompute apparent magnitudes and effective radii. The 
properties of our sample of 48 GCs are summarized in Table lA. 11 
In the Appendix we describe the behaviour of two small er sub- 
samp les, that are obtained by^onsidering either Ferraro et al.l 
d 19991) or lRecio-Blanco et all (120051) . separately. 

3. The "Fundamental Plane" 

After extracting distance moduli \m - M|o and reddening values 
E(B - V) from our sources and reconstructing the apparent ef- 
fective radii and apparent integrated magnitudes m, we can 
proceed to choose the two photometric variables that enter the 
FP from the three definitions (unit-dependent constants, which 
play a role only in determining the zero point, are set by measur- 
ing Re in kpc, in arcsec, and SB,, = 41.40 - 2.5 log(L/2nR^) 
in mag/arcsec^ where L is expressed in units of the solar V-band 
luminositjO): 



log/?, = i|m - M|o+ log -7.31 , 

SBe = OT - 3.1E(B -V)+5 log -H 1.99 , 
M ^m-\m- M|o - 3.1E{B - V) . 



(1) 

(2) 
(3) 



At this stage, some comments on the possible correlations of 
errors are in order We can be quite confident in assuming that 
errors on distance moduli and reddening values do not correlate 
with errors on apparent effective radii and integrated apparent 
magnitudes. On the other hand, as noted earlier, some degree 
of anti-correlation is present between apparent magnitudes and 
effective radii (see Treu et al. 1999, for a discussion of correla- 
tions among analogous parameters of the FP of early type galax- 
ies). From Eq. (|2| above, we see that SB^ is the only photometric 
coordinate of the parameter space that does not depend on dis- 
tance moduli (although it depends on the reddening). Therefore, 
as long as errors on distance moduli dominate the error budget, 
for optimally dealing with the problem of error correlations, it 



- Our source, lTrageretal.Hl995l) . gives surface brightnesses in the V- 
band, although it is actually based on data taken in different passabands. 
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is appropriate to choose from the three photometric parameters 
listed above any of the other two photometric quantities together 
with SBg. Even if the contribution of the error on log r^. to the er- 
ror on SBe is comparable to that of the error on distance moduli, 
it is partly cancelled by the anti-correlated error on the apparent 
magnitude m. 

These arguments confirm that the standard use of the two 
photometric variables \ogRe and SB^ is justified and sensible. 

For an ordinary least square linear regression in the three- 
dimensional parameter space defined by \ogRe, SB^, and log cr, 
we are now left with choosing which coordinate to take as de- 
pendent variable, the other two being regarded as independent 
variables. In the literature the most common choice as depen- 
dent variable is SB^, which is known to prov ide "a better an d 
more stable fit" with respect to other options ( Di org ovskil 1 9951) . 
In fact, different choices give different results, and it is not clear 
which one is the best. This is related to the rather narrow distri- 
bution of the cluster data points on the FP viewed face on (i.e. 
to the datapoints lying more around a line than a plane, see dis- 
cussion in the following sections), a point that we have tested by 
a set of simple simulations. If the data lie approximately on a 
line in parameter space, then any plane passing through that line 
is a good fit to the data. The best fit is different when different 
variables are chosen as dependent variables because the sum of 
squared residuals along that coordinate gets minimized; in this 
way, the fitting algorithm chooses the best-fit plane almost at 
random among the infinitely many planes that pass through the 
line, by fitting the scatter around such line. 

We find the following fit by the Fundamental Plane based on 
SB,: 

SB, = (3.28±0.57)log/?,-(3.59±0.39)logcr-H(27.95+2.94),(4) 
In turn, fitting through log R, we get: 

log/;, = (0.13±0.02)5B,-H(0.30±0.13)logo--(4.94+0.50),(5) 
and, finally, through log cr: 

logo- = -(0.19+0.01)5B,-i-(0.42+0.12)log/;,-H(5.07±0.46).(6) 

In all these fits, we weighted the data-points by the inverse 
squared uncertainty on the relevant dependent variable, neglect- 
ing other uncertainties. The constant term assumes that angular 
effective radii are measured in arcseconds, velocity dispersions 
in kilometers per second, distance moduli and integrated appar- 
ent luminosities are measured in magnitudes, and surface bright- 
ness is expressed in magnitudes per square arcsecond (see also 
the beginning of this Section). 

The relatively large uncertainties on the fit coefficients are 
not the result of overestimating error bars on data points. In fact, 
the reduced for the FP fits is rather high: 9.56, 4.38 and 18.88, 
respectively. These values derive either from underestimating 
observational errors (which is probably not the case, given the 
conservative approach taken in this paper) or, more likely, from 
the presence of intrinsic scatter in the FP relationship. We will 
address this issue further in the next Section, by discussing the 
correlations between the residuals to the FP and other GC obser- 
vational parameters. 

It is customary (see lDiorgovski|[r995h to re-express the FP 
in terms of logR,, even if we prefer, as argued above, to carry 
out the fit through SB,. This involves an algebraic manipulation 
of the coefficients obtained from the SB, based fit. Therefore, 



18 20 
3.28logRo-3.59logcJo + 27.94 



Fig. 2. Edge-on view of the FP for the galactic GCs. 



to get error bars for the derived coefficients it is natural to ap- 
ply standard error propagation techniques, thus obtaining from 
Eq. (H: 

logR, = (1.09+0.31) log cr-H(0.30±0.05)5B,-(8.52±2.35), (7) 

which can be easily compared to the results of other studies 
or to theoretical predictions. We should recall that error bars 
computed through standard error propagation techniques are not 
necessarily accurate when uncertainties are large, because such 
techniques are linear. Figure|2]shows an edge-on view of the FP. 
The scatter in logR, is approximately 0.15. 

Equation & is eas ily compared to previous results, e.g. Eq. 
6 in iDi orgovskil d 1 995h . which we report here in our notation: 

SB, = (2.9 + 0.1)log/?, -(4.1 ±0.2)logo- + (19.8 + 0.1) . (8) 

The coefficients we obtain by fitting through SB, are compat- 
ible, within two sigma, with the results by Diorgovski (19951), 
while the zero point is not, due to different conventions in mea- 
surement units. The results recorded in Eqs. (l4]i-(I7]i are not to 
be seen as interesting because they are consistent with those of 
other investigations, but rather because they demonstrate the dif- 
ficulty in fitting by a plane, thus leading to the discussion of Sect. 

mi 

We finally note that the FP coefficients just found in Eq. (7) 
by fitting through SB, and inverting for log R, are consistent with 
those of the FP f or early-ty pe galaxies, also in the zero-point 
constant, namely (,j0rgensen et al. 1996) : 



logR, = 1.24 log o-o + 033SB, - 8.895 



(9) 



with reported scatter of 0.08 in \ogR„ but the error bars are 
marginally compatible even with the zero-tilt plane: 



log R, = 2.Q log cr + QAQSB, + y . 



(10) 



This result shows that better data are required in order to 
draw convincing physical conclusions. 

3.1. A Fundamental Line? 

If on the FP viewed face-on the data-points are located mainly 
within a thin strip (with width comparable to the scatter around 
the FP) it may be difficult to discriminate between a plane dis- 
tribution and a line distribution. In principle we could argue 
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that, by suitable rescaling of the coordinates, any plane distri- 
bution may be stretched into a line distribution and vice-versa, 
although we should check how the coordinate transformation 
changes the relevant error bars. Rather than following this ap- 
parent paradox as a general problem of statistics, here below we 
adopt the view that the standard choice of coordinates used in 
discussions of the FP has its own physical justification and wish 
to check whether indeed, in the standard FP coordinates, the set 
of points occupies a rather slim cylindrical region of parameter 
space or not. We thus calculate the eigenvalues and eigenvec- 
tors of the tensor of inertia of the unweighted data-points in the 
(SBe - {SBe), log Rf, - (log Rg), log cr - (log cr)) coordinates (i.e. 
in the "center of mass" frame) and obtain for the eigenvalues /i, 
h, and ly. 



h = 1.99 


(11) 


h = 120.68 


(12) 


h = 121.33 


(13) 


and for the respective eij 
dinates): 


jenvectors (in the aforementioned coor- 


VI = (0.98,0.09,-0.15) 


(14) 


V2 = (0.05,0.67,0.74) 


(15) 


V3 = (0.17,-0.73,0.66) 


(16) 



The eigenvector corresponding to the smallest eigenvalue 
sets the direction of the "Fundamental Line". This procedure 
is equivalent to unweighted least square fitting in terms of a 
Fundamental Line, (i.e. finding the straight line which minimizes 
the sum of distances to the points squared) but gives also a mea- 
sure of the degree of axisymmetry of the distribution of data- 
points around such Fundamental Line, in terms of the difference 
between the other two eigenvalues h and ly. 



R = 



2\h - h\ 
(h + h) 



= 0.005, 



(17) 



which is remarkably small. 

Indeed, it was already noted (iBurstein et al.ll 19971 iBellazziriil 

119981) that the galactic GCs lie close to a straight Une in the 
log Re, SBe, logo- parameter space. But, to the best of our 
knowledge, the present paper is the first to identify and quan- 
tify the degree of axisymmetry around such line, a fact that 
demands an adequate theoretical explanation. Since IBellazziriil 
( Il998l) suggests that "at earlier times, globular clusters popu- 
lated a line in the three-dimensional S-space, i.e. their origi- 
nal dynamical structure was fully determined by a single phys- 
ical parameter", we also looked for coiTelations between resid- 
uals to the Fundament al Line and GC relative age (taken from 
iDe Angeh et al.l l2005h . finding only a very small, but posi- 
tive, correlation coeflicient of 0.05, which further decreases to 
0.03 if more robust non-parametric correlation estimators (such 
as the Spearman estimator (e.g., see IWassermanI 120061) ) are 
used. However, this apparently negative result is not conclusive, 
because our sample contains mainl y metal-poor clusters (see 
Sect. |5j, and 'De Ange h et all (l2005h clearly state that "the age 
dispersion for the metal-poor clusters is 0.6 Gyr (rms), consis- 
tent with a null age dispersion". In conclusion, to clarify this 
point a larger sample would be desired. 



3.2. A pure photometric scaling law? 

We note that a Fundamental Line relation of the type identified 
above would imply the existence of a pure photometric scaling 
law approximately of the form 



log Re -^^^e 



LocR 



-2 



(18) 



(19) 



This result is quite unexpected because the trend suggeste d by 
Eg. (fT9T l would be qualitatively opposite to that of the iFishI 
('1964') law. In addition, we should recall that, on a sample of 
143 GCs, Diorgovski & Mevlan ( 1994) conclude that no signif- 
icant correlation exists between integrated absolute magnitude 
and the logarithm of the effective radius. From the relevant plot 
in their Fig. 6, we note that an underlying trend might be masked 
by the presence of two or three outlier points; furthermore, in 
their Table 1 Diorgovski & Mevlan ( 1994) actually report a pos- 
itive correlation (with r - 0.819) between logRe and SBe ((yUy);, 
in their notation), which can be recast into a scaling law of the 
form of Eq. (fT9] l. recovering the correct sign for the exponent. 

On the other hand, the correlation that we find in the pho- 
tometric plane turns out to be in qualitative agre ement with 
the ge neral trer id noted for elhptical galaxi e s by iKormendv 



tie ge neral trer ia notea tor elliptical galaxi e s dv lJs.ormena\ 
(11977 ): see also lHamabe & Kormendvl(ll987h : ICapaccioli etal 
(Il992h . We have checked that, for the sa mple of 



GCs ob- 



tained b y intersecting the tw o samples of Ferraro et al. (1999| 
and Recio-Blanco et al] (12005) with that of Trager et al. (19951 
the correlation that we find in the photometric plane indeed has a 
consistent slope but a diff'erent zero point with respect to the rela- 
tion Hev - 2.94logRe + 19.48 reported by Hamabe & Kormendy. 
To be sure, the fiev variable is very sensitive to the quality of 
the surface brightness profiles; therefore, a direct deta i led co m- 
parison with the least squares regression of lKormendvl ( Il977h in 
the ij.ev variable cannot be performed easily and would require 
additional discussion. 

In fact, the value of the exponent in Eq. ( fT9] ) is very sensi- 
tive to the adopted fitting procedure. The projected Fundamental 
Line is not equal to the straight line that minimizes the sum of the 
squared residuals along one coordinate in the plane (log Re, SBe), 
i.e. it does not coincide with the line one would get with simple 
linear least squares fitting. To recover the same relation as pre- 
dicted by the projected Fundamental Line, one should perform 
an impartial regression in the (log Re, SBe) plane, i.e. one should 
minimize the sum of squared distances to the line, not just the 
residuals in one coordinate. Moreover, it should be noted that the 
least squares method implicitly used for fitting the Fundamental 
Line is non-robust, i.e. its results may be excessively influenced 
by outliers; this is obviously true also for the Fundamental Line 
projection and the derived scaling relation. 

The scatter around the Fundamental Line has the poten- 
tial for significantly blurring the univariate correlation ob- 
tained by projecting the Fundamental Line on the photometric 
{log Re, SBe) plane. Therefore, we anticipate that either a nar- 
rowing of the sample to a smaller set of GCs or the identification 
of a variable 77 capable of reducing that scatter (77 could be the 
central slope of the surface bri ghtness profile, as in the case of 
the FP, or age as suggested b v Bellazzinil 11998 *) might lead to 
a much improved univariate relation between log Re and SBe. 
Possibly, a new bivariate correlation could be sought for in the 
[logRe, SBe, 77) space. 



Pasquato & Berlin: Fundamental Plane of Globular Clusters 



7 



A clear interest in a scaling law of the form of Eq. ( fT9l ) lies 
in the fact that its study does not require the acquisition of kine- 
matical information. Therefore, we should be able to test it on 
wider samples and possibly on globular cluster systems in exter- 
nal galaxies. This will soon be addressed in a follow-up paper 
The exclusion of the log cr coordinate would also allow us to 
consider samples definitely more homogeneous from the statis- 
tical point of view. Finally, the existence of a pure photometric 
scaling relation suggests that the physical basis for the empirical 
scaling should be sought more in the context of stellar popula- 
tions than within the virial theorem and dynamical arguments. 

4. Two unexpected features 

4. 1 . Correlation of residuals with the central slope of the 
surface brightness profile 

The scatter around the FP, as quantified by the value of the re- 
duced x^, significantly greater than unity, might be interpreted 
as an indication of the role of nonlinearities or of a fourth pa- 
rameter. 

First of all, we have looked for residual correlations with 
the coordinates defining the space in which the FP lives, namely 
log Rf, SBg, and log cr. By looking at the relevant scatter plot, the 
residuals show no recognizable pattern with either log R^, or SBe, 
or log cr; from a qualitative point of view, this suggests that non- 
linear corrections to the FP, if present, are unobservable, at least 
with present-day data quality. Quantitatively, the null hypothe- 
sis that the coefficients of nonlinear terms in the form (log/J^,)^, 
log Rf, log cr, . . . , be equal to cannot be refuted even at a 20% 
significance level. 

Then we have moved on to consider correlations with other 
quantities, mainly of dynamical interest. In particular, we have 
referred to: King model concentration parameter, central slope 
of the surface brightness profiles (as defined and measured by 
iNovola & Gebhardt,2006.) . position in the Galaxy (i.e. galacto- 
centric distance and vertical distance from the galactic disk), 
age, relaxation time referred to the effective radius, metallic- 
ity [Fe/H], and colour B - V. In this respect, our investi- 
gation is not complete, because other quantities, in particular 
colour-magnitude-diagram related quantities, might have been 
addressed. 

Consisten t with indicatio ns from previous studies (e.g., see 

iDiorgovski & MevlanI 1 19941) which excluded correlations be- 
tween GC properties, even in the case of FP residuals no sig- 
nificant correlations are found, with just one interesting excep- 
tion, that is a statistically significant correlation of the FP resid- 
uals with the central slope of the surface brightness profile. In 
view of the results obtained by McLaughlin (2000) (see also 
iBarmbv et alj|2007l) . a trend with galactocentric position would 
be expected; the fact that the present analysis is unable to detect 
it could be due to the selection effect in galactocentric distance 
characterizing our sample. Figure [3] shows a plot of FP residual s 
against the central slopes taken from Novola & Gebhardtl(l2006h . 
The trend of smaller effective radii with increasing slopes is ev- 
ident from the scatter plot; the value of the related correlation 
coefficient is 0.73. The presence of such a correlation reinforces 
the view that the origin of the FP scatter is indeed intrinsic. 
This result is particularly interesting, because it relates central 
to global prope rties of these stellar systems, much like the Mbh- 
o- relationship dFerrarese & Merrittll2000l: iGebhardt et al.ll2000l) 
does for bulges and elliptical galaxies. Indeed, an analogous ar- 
gument has been proposed by Faber et akl ( Il997h in their study 
of the core regions of elliptical galaxies and spiral bulges; in 
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Fig. 3. Correlation of the FP residuals with slope of the cen- 
tral surface br ightness profile as defined and measured by 
iNovola & Gebhardt (2006.) . 



particular, iFaber et akl (Il997h state that "Cores follow a funda- 
mental plane that parallels the global fundamental plane for hot 
galaxies but is 30% thicker Some of this extra thickness may be 
due to the effect of massive black holes (BHs) on central velocity 
dispersions". We note that in some GCs the gravitational sphere 
of influence of an intermediate mass black hole ( f or exa mple, 
such a black hole has been argued bv lNovola et al.l (l2008l) to be 
present in Omega Cen) could have an angular size comparable to 
the slit size used for integrated light measurements of the veloc- 
ity dispersion. This effect might be at the origin of the observed 
correlation; fitting cuspy GC cores with double-power law pro- 
files, (see Byunet al. 1996) or with simple dynamical models 
(Ci pollina & Bertinlll994t) might shed further light on this ques- 
tion. 

We have also checked whether the scatter in the FP can be 
related to the presence, in our sample, of objects that are sus- 
pected t o be non-genuine galactic GCs, su ch as members of GC 
streams ( Lvnden-Bell & Lvnden-B ell 1995), which are probable 
remnants of GC systems tidally stripped from a smaller galaxy 
accreted by the Milky Way. We have compared the distribution 
of residuals t o the FP for those GCs that are part of such streams 
(according to ! Gao et al.l2007l) with that of the rest of our sample 
through a Kolmogorov-Smirnov test and found no evidence of 
significant differences. Therefore, we argue that GC streams are 
not responsible for the excessive FP scatter, or, in other words, 
that, in this context, stream GCs are dynamically similar to na- 
tive GCs. 



4.2. Distribution of effective virial coefficients 

We have reconstructed the distribution of Kv/{M/L) = 
GL/(ReCr^) (virial coefficient divid ed by the mass-to- light ra- 
tio; for a thorough discussion, see iBertin etalJl2002h for our 
sample, making use of kernel density estimation techniques (see 
Wasserman 2006). In many respects, these non-parametric tech- 
niques are an extension of histograms, best suited to spot struc- 
ture that histograms may hide. The Kv/{M/L) distribution looks 
bimodal, even if we choose a slightly larger bandwidth than op- 
timal (1.2 times) so to err on the side of oversmoothing (under- 
smoothing could easily produce fictitious bi- or multi-modal dis- 
tributions). The indication of bimodality is stronger if the anal- 
ysis is restricted to the 28 systems for which the line-of-sight 
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Fig. 4. Probability density distribution of log Kyi {Ml L) for clus- 
ters with velocity dispersions measured in integrated light (solid 
line) and for the entire sample (dashed line), reconstructed by ap- 
plication of the non-parametric method of kernel density estima- 
tion. Units have been chosen in such a way that Kyi {MIL) - 1 
forNGC 104. 



velocity dispersion is obtained from integrated light measure- 
ments, while for the entire sample of 48 GCs it almost disap- 
pears. Figure ID illustrates the reconstructed probability density 
distribution of log Kyi {MIL) for the restricted and the unre- 
stricted case, with the restricted case exihibiting the superpo- 
sition of a narrow peak around a negative value onto a broader 
distribution with positive mean. Therefore we note that inhomo- 
geneity of velocity dispersion measurements is not responsible 
for the bimodality, which appears more clearly when the sam- 
ple is restricted to integrated light velocity dispersion measure- 
ments only. This result is also suggestive of an interpretation in 
terms of the presence of intermediate mass black holes, which 
would naturally affect the value of the central velocity disper- 
sion (best diagnosed by the integrated light measurements), in 
some of the globular clusters, but, as for the feature noted in the 
previous subsection, a thorough investigation on a larger sam- 
ple would be required in order to make a fully convincing case. 
Even though the mass-to-light ratio of a stellar population may 
be heavily influenced by its metallicity, we feel confident in ex- 
cluding a metallicity-dependent effect as responsible for the bi- 
modality of the log Kvl {MIL) because no correlation is observed 
between eff'ective virial coefficient and metallicity; using [FelH] 
data from the McMaster Catalogue, we find a correlation coeffi- 
cient of 0.04 between these two quantities, which is statistically 
compatible with zero. 

In addition, this result makes a point which is clearly against 
the zero tilt FP scenario, by showing that the log Kyi {MIL) val- 
ues have a non-trivial distribution, i.e. inconsistent with a sharply 
peaked distribution broadened into a bell-shaped curve by obser- 
vational errors. 



5. Selection effects 

As noted earlier, the 
Globular Cluster Parameters 



McMaster Catalo; 
jHarrisI \l99' 



ue of Milky Way 
lists 151 globular 
clusters, while the sample considered in the present paper is 
much s maller, comprising 48 GCs. By comparison, the sample 
used by Diorgovskil d 19951) contains 55 globular clusters, those 
for which a measurement of line-of-sight velocity dispersions is 



available in iPrvor & MevlanI ( Il993h . Under these conditions a 
thorough study of selection effects is mandatory. 

We proceeded as follows. First of all we defined a list of rel- 
evant parameters, extracted from those tabulated in the 2003 ver- 
sion of the McMaster Catalogue, namely the distance of the GC 
from the Galactic Center Rgc, the vertical distance to the Galactic 
Plane Z, the distance from the Sun do, the absolute integrated 
magnitude My, the integrated colour B-V (uncorrected for red- 
dening), the metallicity as measured by the [FelH] parameter, 
the King model concentration parameter c, the projected effec- 
tive radius R^, the logarithm of the relaxation time referred to the 
effective radius log fg, and the central surface brightness fiy. We 
then compared the distribution of these parameters for our sam- 
ple with that of the complementary sample, that is of the sam- 
ple containing all of the GCs in the McMaster catalogue except 
those of our sample. 

In assessing the importance of selection effects, we used the 
quantities taken from the McMaster Catalogue; we did not try to 
recompute any of them anew for consistency reasons, because 
we would otherwise treat our sample (for which quality data are 
available) in a different way with respect to the rest of the clus- 
ters. 

The comparison has been performed through the non- 
parametric Kolmogorov-Smirnov test. The test gives an estimate 
of the probability of the sample and of its complementary set 
being drawn from the same distribution. If this probability is ex- 
tremely low, say less than 10"^, we can quite confidently assert 
that a selection effect is present, i.e. that the sample we consid- 
ered is not statistically representative of the whole population of 
galactic GCs. Being non-parametric, this test does not assume 
a given functional dependence for the distribution function of 
the data, which is unknown in principle; this is an important re- 
quirement in this context, because the distribution functions of 
the various parameters we considered are often far from gaussian 
and may be highly skewed or have long tails. 

We performed a similar test also for the wider sample of 
55 GC s with veloci t y disp ersions measured by iPrvor & MevlanI 
(ll993h . lDiorgovskil (11995') assumes no relevant selection effects 
on that sample, at variance with our analysis, that shows that 
some selection effects are indeed present at a significant level. 

Our statistical tests are summarized in Ta b. 2 (our sample 
against the rest of galactic GCs) and Tab.[3](the lPrvor & MevlanI 
(1993) sample against the rest of galactic GCs), which also list 
the mean and standard deviation of the chosen sample and of its 
complementary sample. 

The results of our statistical analysis can be summarized as 
follows. Selection effects are definitely present, although they 
do not necessa rily hinder the con clusions of the present paper 
Both samples (iDiorgovskil (Il995l) and ours) are biased in favor 
of clusters relatively close to the Sun; both are biased towards in- 
trinsically more luminous and centrally brighter clusters. These 
selection effects, although expected, must be kept in mind when 
addressing issues such as whether the FP properties change with 
Galactocentric distance. 

On the other hand, it is somewhat surprising to find that 
both samples are biased in the direction of bluer and more cen- 
trally concentrated clusters (even though this latter feature may 
partly result from the fact that the McMaster Catalogue assigns 
a concentration parameter of 2.5 to GCs recognized to be in 
the Post Core Collapse phase of their evolution). To be sure, 
we are aware that concentration parameters show a trend with 
GC integrated luminosities (Diorgovski & Meylan 1994, see), 
so the bias towards concentrated clusters might be a byproduct 
of the selection effect in magnitude. Our sample has somewhat 
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Table 2. Comparison of the 48 GC sample of this paper with the 
complementary sample within the McMaster Catalogue, from 
which the data listed in this Table are taken. Kolmogorov- 
Smirnov p-values are approximated to one significant digit, 
while all other numbers are truncated to the second digit past 
the decimal dot. 



ID 


Pk&s " 








0-/ 


R^, (kpc) 


0.2 


10.05 


12.96 


13.62 


22.02 


Z (kpc) 


0.1 


1.94 


0.53 


8.13 


18.09 


do (kpc) 


0.2 


11.72 


16.80 


12.41 


20.60 


Mv 


0.0002 


-7.78 


-6.69 


1.07 


1.58 


B-V 


0.004 


0.90 


1.18 


0.27 


0.48 


[Fe/H] 


0.1 


-1.45 


-1.22 


0.49 


0.58 


c 


0.0006 


1.82 


1.44 


0.54 


0.54 


Re (pc) 


0.04 


3.63 


4.94 


3.22 


4.68 


log h (yr) 


0.2 


9.09 


9.04 


0.39 


0.54 




8 • 10-' 


16.99 


19.82 


2.01 


2.66 



" p- value for the null hypothesis of the sample and its complementary 

sample coming from the same distribution 

* Mean value of the relevant parameter on the sample 

Mean value of the relevant parameter on the complementary sample 

RMS of the relevant parameter on the sample 
'' RMS of the relevant parameter on the complementary sample 



Table 3. Comparison of the 55 GC sample of Djorgovski 
(1995) with the complementary sample within the McMaster 
Catalogue, from which the data listed in this Table are taken. 
Kolmogorov-Smirnov p-values are approximated to one signif- 
icant digit, while all other numbers are truncated to the second 
digit past the decimal dot. See Table 2 for column heading ex- 
planations. 



ID 


Pk&s 










Rgc (kpc) 


0.1 


9.70 


13.38 


12.90 


22.72 


Z (kpc) 


0.5 


1.54 


0.66 


7.71 


18.73 


do (kpc) 


0.1 


11.33 


17.39 


11.94 


21.13 


Mv 


2 • 10"^ 


-7.85 


-6.57 


1.15 


1.52 


B-V 


0.002 


0.90 


1.21 


0.26 


0.49 


[Fe/H] 


0.02 


-1.48 


-1.19 


0.46 


0.59 


c 


0.0009 


1.79 


1.43 


0.52 


0.56 


Re (pc) 


0.04 


3.62 


5.05 


3.04 


4.84 


log/,, (yr) 


0.2 


9.11 


9.02 


0.40 


0.54 




2 • lO"** 


16.96 


20.06 


1.94 


2.58 



lower mean metalUcity than average; this selection effect be- 
com es even stronger if on ly the subsample with distance moduli 
from lFerraro et al.l(ll999h is used (se e Appendix). 

The clusters in the iDiorgovskil (fl995) sample are slightly 
smaller than average, in terms of while for both samples 
there seems to be no bias in the relaxation time referred to Re- 
Therefore, we feel confident in ex cluding major dyna mical dif- 
ferences between the two samples ('Diorg ovskil (1 19951) and ours) 
and the overall population of galactic GCs, which makes it rea- 
sonable to extend our conclusions to the whole system of galac- 
tic GCs. 

In closing this Section we should mention that selection ef- 
fects might be present in more than one variable taken together, 
even when they do not show up by considering only one variable 
at a time, as we did. That is, if x and y are two variables, the dis- 
tribution function of our sample could be diff'erent from that of 



the rest of the GCs on the (x, y) plane, while being indistinguish- 
able when projected onto the x or the y axes. This issue has been 
considered to be beyond the scope of the present paper. 

6. Summary, discussion, and conclusions 

6.1. Results 

Based on the identification of an optimal sample of 48 globular 
clusters selected in terms of homogeneity criteria for the photo- 
metric data available in the literature, the most important results 
of the present study are the following: 

- We have determined the coefficients of the FP for the galactic 
globular cluster system, with error bars (see Eq.|7]). 

- Since the reduced of the fit has been found to be signif- 
icantly greater than unity, we conclude that either we have 
underestimated eiTors or, more likely, that the relatively large 
scatter around the FP (of about 0.15 in logR,,) is intrinsic. 
We found that the so-called GC streams are likely to be not 
responsible for the scatter around the FP. 

- We have shown that in the standard FP coordinates the set 
of points representing our sample of GCs occupies a rather 
slim cylindrical region of parameter space, which suggests 
that the relevant scaling relation might be around a line, 
rather than a plane; this confirms a result noted previously by 
Bellazzini (1998). We think that this is the origin of the dif- 
ficult i es in fits by a plan e noted earlier (e.g., see lDjorgovskil 
[T99llDe MicheUllIOOSi) . In addition, we have derived eigen- 
values and eigenvectors for the tensor of inertia of the dis- 
tribution of data-points in the SBg, logR,,, log cr coordinates 
and found that such distribution is remarkably axisymmetric 
around a line, mainly along SBf,. 

- For our sample, we found no statistically significant corre- 
lation of FP residuals with metallicity, colour, concentration 
parameter, position in the Galaxy, ages, and relaxation times, 
but we found a potentially interesting correlation with the 
central slope of the surface brightness profile. 

- By means of kernel density estimation techniques, we have 
reconstructed the probability density distribution of the ef- 
fective virial coefficient Kv/{M/L) of GCs, providing evi- 
dence for bimodality. We argue that this feature is not related 
to the standard disk/halo dichotomy, because our sample 
comprises only clusters with metallicity below about -0.8, 
which are therefore likely to be part of the halo population. 

- By means of non-parametric statistical tests, we studied the 
selection effects aff ecting our sample and the sample used by 
IDiorgovskil (1 19951) in his pioneering paper on the FP for the 
galactic GCs, demonstrating that some selection effects are 
present (in particular, our sample is biased in the direction of 
bright, metal poor clusters close to the Sun). No significant 
selection effects are present in the relaxation time refeiTed 
to the effective radius, which is reassuring because it sug- 
gests that our sample is not dynamically different from the 
remaining population of galactic GCs. 

6.2. Discussion and conclusions 

The eiTor bars on the FP coefficients are rather large. Therefore, 
on the basis of currently available data it is difficult to tell 
whether indeed the FP relation of GCs is a continuation of that 
of early-type galaxies. The limited luminosity range of galactic 
GCs may be an issue here; studies of extragalactic GC systems 
such as that of NG C 5128, which reach a higher upper lum inos- 
ity hmit (e.g., see IReikuba et al.ll2007l; [Martini & Holl2004 . are 
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in line with the present study. In principle, the current FP fit does 
not exclude (at least at a marginal level) that a no-tilt relation is 
at the basis of the observed correlations. On the other hand, we 
found concrete indications that the shape of the distribution of 
Ky/iM/L) is bimodal; if such indications will be confirmed, this 
would clearly point against the naive interpretation that the FP 
just reflects the virial constraint. 

The hints of a bimodal distribution of the effective virial co- 
efficient Kv/(M/L) brought us to consider one specific physi- 
cal ingredient that might play an important role in the statistical 
properties of the GC system, that is the possible presence of in- 
termediate mass black holes. In fact, such presence would also 
have a counterpart in the second unexpected feature noted in our 
statistical investigation, i.e., the trend of FP residuals with the 
central slope of surface brightness profiles. 

As to the possibility that the Fundamental Plane should be 
replaced, for globular clusters, by a Fundamental Line relation, 
before trying to put forward a physical interpretation we prefer 
to wait for confirmations from a wider sample or from studies 
of globular cluster systems in external galaxies. The possibility 
of getting soon an answer to this issue is reasonable because, 
as discussed at the end of Sect. 3, the Fundamental Line relation 
appears to imply the existence of a pure photometric scaling law. 
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Appendix A: Discussion of tlie FP on two 
subsamples witli uniform distance moduli 

For the sub- s ample of 35 GCs with distance moduli from 
iFerraro et al.l d 19991) we find the following fit by the 
Fundamental Plane based on SB^: 

SB, = (3.2 + 0.7) log 7?, -(3.4 + 0.5) logo- -I- (27.96 ±3.5) .(A.l) 
In turn, fitting through log we get 

\ogRe = (0.14±0.03)5B,+(0.3±0.1)logo--(5.04 + 0.6),(A.2) 
and, finally, through log a". 

logo- = -(0.18±0.02)5B, + (0.4±0.2)log7;, + (5.15±0.6).(A.3) 

In tu rn, for the sub-sarnple of 34 GCs with distance moduli 
froml ecio-Blanco et al.l j2005h we find the following fit by the 
Fundamental Plane based on SB/. 

SB, = (2.66+0.75) log7;,-(4.02±0.47) log o-+(31.70+3.9).(A.4) 
In turn, fitting through log R, we get 

log/?, = (0.12+0.03)5fi,+(0.26+0.18) log o--(5. 19+0.65), (A.5) 
and, finally, through log cr: 



These results are compatible with those obtained using the 
global sa mple of 48 GCs co nsidered in the main text of the pa- 
per Since iRecio-Blanco et al. (2005) and Ferraro et al. (1999!) 
use a different definition of the ZAHB level, these results show 
that our conclusions (and in particular the FP coefficients) are 
relatively independent of the details of ZAHB fitting. 

We have also checked that the general results listed in Sect. 
l6.1l hold for the two subsamples if studied separately. 
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ID 






log r/ 


\M - m\l 


6\M- m\l 


E{B - 


t/logE/c/log 


Distance sample'' 


NGC 104 


9.72 


4.00 


2.22 


13.27 


0.08 


0.04 


0.00 ± 0.04 


R05 F99 


NGC 288 


2.61 


8.45 


2.16 


14.67 


0.30 


0.03 




F99 


NGC 362 


6.15 


6.26 


1.69 


14.69 


0.10 


0.05 




R05 F99 


NGC 1851 


10.21 


7.32 


1.46 


15.44 


0.09 


0.02 


-0.38 ±0.11 


R05 F99 


NGC 1904 


4.59 


7.70 


1.93 


15.67 


0.08 


0.01 


-0.03 ± 0.07 


R05 F99 


NGC 2419 


2.66 


10.46 


1.17 


19.88 


0.29 


0.03 




F99 


NGC 2808 


14.04 


6.42 


1.77 


14.96 


0.13 


0.23 


-0.06 ± 0.07 


R05 F99 


NGC 3201 


4.30 


7.21 


2.30 


13.58 


0.13 


0.21 




R05 F99 


NGC 4147 


2.54 


10.32 


1.46 


16.42 


0.10 


0.02 




R05 F99 


NGC 4590 


2.39 


8.32 


1.96 


15.16 


0.07 


0.04 




R05 F99 


NGC 5053 


1.15 


9.79 


2.29 


16.17 


0.19 


0.03 




F99 


NGC 5272 


5.22 


6.65 


1.83 


14.99 


0.14 


0.01 


-0.05 + 0.10 


F99 


NGC 5286 


8.00 


7.41 


1.90 


15.25 


0.29 


0.24 


-0.28 + 0.11 


F99 


NGC 5466 


1.57 


9.11 


2.13 


16.12 


0.29 


0.00 




F99 


NGC 5694 


5.50 


10.17 


1.43 


17.80 


0.11 


0.09 


-0.19 + 0.11 


R05 F99 


NGC 5824 


10.73 


8.90 


1.44 


17.63 


0.11 


0.13 


-0.36 + 0.16 


R05 F99 


NGC 5904 


5.24 


5.93 


1.97 


14.44 


0.09 


0.03 


0.05 + 0.07 


R05 F99 


NGC 5946 


3.70 


9.74 


1.62 


15.41 


0.29 


0.55 




R05 


NGC 6093 


13.32 


7.47 


1.56 


15.09 


0.15 


0.18 


-0.16 + 0.07 


R05 F99 


NGC 6121 


3.59 


5.67 


2.34 


11.62 


0.30 


0.36 




F99 


NGC 6171 


3.67 


8.29 


2.16 


14.07 


0.17 


0.33 




R05 F99 


NGC 6205 


6.62 


6.00 


1.96 


14.45 


0.08 


0.02 


-0.10 + 0.15 


R05 F99 


NGC 6218 


3.92 


7.27 


2.02 


13.66 


0.23 


0.18 




R05 F99 


NGC 6254 


5.48 


6.70 


2.08 


13.35 


0.29 


0.28 


0.05 + 0.07 


F99 


NGC 6256 


6.50 


10.91 


1.89 


14.97 


0.37 


0.84 




R05 


NGC 6266 


14.23 


6.50 


1.81 


14.31 


0.21 


0.47 


-0.13 + 0.08 


R05 F99 


NGC 6284 


6.20 


8.82 


1.70 


16.08 


0.21 


0.28 


-0.55 + 0.14 


R05 


NGC 6293 


7.60 


8.28 


1.76 


14.88 


0.24 


0.39 


-0.67 + 0.08 


R05 


NGC 6325 


5.80 


10.30 


1.88 


14.82 


0.39 


0.89 




R05 


NGC 6341 


5.00 


6.53 


1.69 


14.74 


0.29 


0.02 


-0.01 + 0.04 


F99 


NGC 6342 


4.60 


9.80 


1.68 


15.00 


0.26 


0.44 




R05 


NGC 6366 


0.93 


9.27 


2.26 


12.88 


0.30 


0.69 




F99 


NGC 6362 


2.76 


7.71 


2.24 


14.45 


0.16 


0.09 




R05 


NGC 6388 


18.90 


6.95 


1.54 


15.38 


0.24 


0.38 


-0.13 + 0.07 


R05 


NGC 6397 


3.48 


6.06 


2.28 


11.89 


0.29 


0.18 


-0.37 + 0.11 


F99 


NGC 6441 


17.98 


7.35 


1.50 


15.82 


0.26 


0.45 


-0.02 + 0.12 


R05 


NGC 6522 


6.70 


8.34 


1.71 


14.77 


0.28 


0.50 




R05 


NGC 6535 


2.15 


10.53 


1.66 


14.29 


0.45 


0.32 


-0.50 + 0.18 


F99 


NGC 6624 


5.54 


8.00 


1.66 


14.45 


0.21 


0.27 


-0.32 + 0.16 


R05 


NGC 6681 


8.21 


8.12 


1.68 


14.99 


0.10 


0.07 


-0.82 + 0.09 


R05 F99 


NGC 6712 


4.03 


8.03 


2.03 


14.15 


0.15 


0.46 


0.02 + 0.05 


R05 F99 


NGC 6752 


4.50 


5.75 


2.06 


13.13 


0.45 


0.04 


-0.03 + 0.15 


F99 


NGC 6809 


4.12 


6.96 


2.21 


13.78 


0.29 


0.07 




F99 


NGC 6864 


10.30 


8.51 


1.46 


16.63 


0.18 


0.16 




R05 


NGC 6934 


4.91 


8.86 


1.62 


16.03 


0.11 


0.11 




R05 F99 


NGC 7078 


12.98 


6.42 


1.78 


15.19 


0.11 


0.09 


-0.66 + 0.11 


R05 F99 


NGC 7089 


7.39 


6.46 


1.85 


15.35 


0.15 


0.05 


0.05+0.11 


R05 


NGC 7099 


5.13 


7.48 


1.93 


14.72 


0.08 


0.03 


-0.57 + 0.11 


R05 F99 



° km/s. lPryor & Meylanl (Il993h 
V-band mag, this paper 

arcse c, this paper 

mag.lFerraro et al.l (IT9 99) and'Recio-Blanco et al.' ('2005") 
" iFerraro et al. I (119991): iRecio-Blanco et al.. (,200 5) and this paper (see text) 
^ mag. lRecio-Blanco et al.1 (l2005h and iHarrisI (11996) 
^ dimensionless,'Noyola & Gebhardt (2006^ 

F99 indicates that the GC is listed bv .Feiraro et al.l(ll999l) . R05 that it is listed bv lRecio-Blanco et al.1 (l2005l) . 



